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Abstract 

We study Mach shocks generated by fast partonic jets propagating through a deconfined 
strongly-interacting matter. Our main goal is to take into account different types of collective 
motion during the formation and evolution of this matter. We predict a significant deformation 
of Mach shocks in central Au+Au collisions at RHIC and LHC energies as compared to the case 
of jet propagation in a static medium. The observed broadening of the near-side two-particle 
correlations in pseudorapidity space is explained by the Bjorken-like longitudinal expansion. 
Three-particle correlation measurements are proposed for a more detailed study of the Mach 
shock waves. 

PACS numbers: 25.75.-q, 25.75.Ld, 47.40.-x 
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I. INTRODUCTION 

Sideward peaks have been recently observed in azimuthal distributions 

of secondaries associated with the high hadrons in central Au+Au collisions at 
y/ $nn = 200 GeV. In Ref. Q] such peaks were predicted as a signature of Mach shocks 
created by partonic jets propagating through a quark-gluon plasma (QGP) formed in 



a heavy-ion collision. Ana 



hadronic matter 



ogous Mach shock waves were studied previously in a cold 



10l | as well as in nuclear Fermi liquids 
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12j. The Mach- 

shock induced electron emission from metal surfaces have been predicted lldl and then 

n 

observed experimentally in Ref. [1J|. Recently, Mach shocks from jets in the QGP have 
been studied in Ref. Il_5|| by using a linearized fluid-dynamical approach. It has been 

nn 

argued in Refs. p, H6J that Mach-like motions of quark-gluon matter can appear via the 
excitation of collective plasmon waves by the moving color charge associated with the 
leading jet. 

Most of these studies were dealing with the idealized case of homogeneous, static mat- 
ter. On the other hand, collective expansion effects are known to be strong in relativistic 
colons of heav, ouolei 0. Fo r ex^nple, thecal Us of RH.C data give for most central 
events the radial flow velocities Vf ~ 0.6c. The flow effects accompanying jet propaga- 
tion through the expanding QGP have been considered within different approaches in 
Refs. Q,QQ- 

Numerical solutions of fluid-dynamical equations with an additional 
external source have been studied in Ref. |2l| . 

It is well known [22] that a point-like perturbation (a small body, a hadron or parton 
etc.) moving with a supersonic speed in the spatially homogeneous ideal fluid produces 
the so-called Mach region of the perturbed matter. In the fluid rest frame (FRF) the 
Mach region has a conical shape with an opening angle with respect to the direction of 
particle propagation given by the expression j^] 



'M 

" V 



where c s denotes the sound velocity of the unperturbed (upstream) fluid and v is the 
particle velocity with respect to the fluid. In the FRF, trajectories of fluid elements 
(perpendicular to the surface of the Mach cone) are inclined at the angle A9 = n/2 — 9 M 



with respect to v . Strictly speaking, formula is applicable only for weak, sound-like 
perturbations. It is certainly not valid for space-time regions close to a leading particle. 
Nevertheless, we shall use this simple expression for a qualitative analysis of flow effects. 
Following Refs. j^, Q] one can estimate the angle of preferential emission of secondaries 
associated with a fast jet in the QGP. Substituting 28] v = 1, c s = l/y/3 into Eq. (0) 
gives the value A8 ~ 0.96 . This agrees well with positions of maxima of the away-side 
two-particle distributions observed in central Au+Au collisions at RHIC energies. 

In this paper we study influence of flow effects on the properties of Mach shocks cre- 
ated by a high-energy parton in expanding QGP. In Sect. II we apply simple kinematic 
considerations to study weak Mach shocks propagating in transverse and collinear direc- 
tion with respect to the flow velocity. In Sect. Ill we use the method of characteristics to 
calculate the width of the perturbed region in the rapidity space. Finally, in Sect. IV we 
summarize our results and give outlook for future studies. 



II. DEFORMATION OF MACH SHOCKS DUE TO RADIAL FLOW 

Let us consider first the case when the away-side jet propagates with velocity v parallel 
to the matter flow velocity u . Assuming that u does not change with space and time, 
and performing the Lorentz boost to the FRF, one can see that a weak Mach shock has 
a conical shape with the axis along v . In this reference frame, the shock front angle 9m 
is given by Eq. (0). Transformation from the FRF to the cm. frame (CMF) shows that 
the Mach region remains conical, but the Mach angle becomes smaller in the CMF: 

tan 6 M = — tan 6 M , (2) 

where r y u = (1 — w 2 ) -1 / 2 is the Lorentz factor corresponding to the flow velocity u. 
Eqs. (£Q)-(J21) give the resulting expression for the Mach angle in the CMF 

(3) 




where 

~ v =F u 



1 =F vu 
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FIG. 1: Mach cone angles for jet propagating collinearly to the matter flow as a function of fluid 
velocity u . Different curves correspond to different values of sound velocity c s . 

and upper (lower) sign corresponds to the jet's motion in (or opposite to) the direction 
of collective flow. For ultrarelativistic jets (v — > 1) one can take v ~ 1 which leads to a 
simpler expression 

(, ^ tan " l ( f ?) = sin " 1 ( c -/i^)' (5) 

where j s — (1 — c 2 ) -1 / 2 . According to Eq. in the ultrarelativistic limit 9m does not 
depend on the direction of flow with respect to the jet. The Mach cone becomes more 
narrow as compared to jet propagation in static matter. This narrowing effect has a purely 
relativistic origin. Indeed, the difference between 9 m from Eq. (jSJ) and the Mach angle 
in absence of flow (lim 9m = sin -1 c s ) is of the second order in the collective velocity u . 
The Mach angle calculated from Eq . (15)) is shown in Fig. ^ as a function of u for different 
sound velocities c, . Following Ref. H, the value ef = 1/5 is identified with the hadrome 
matter and c 2 s = 1/3 with ideal QGP composed of massless quarks and gluons. The value 
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= 2/3 is chosen to represent a strongly coupled QGP 

The case of a jet propagating at nonzero angle with respect to the flow velocity is more 
complicated. As will be shown later, Mach shocks become nonconical for non-collinear 
flows. For simplicity, below we study only the case when the jet and flow velocities are 
orthogonal to each other, v J_ u . Let axes OX and OY be directed along u and v , 
respectively. As before, we first make transition to the FRF by performing the Lorentz 
boost along the OX axis. The components of the jet velocity v in the new reference frame 
are equal to 



)yl — it 2 . (6) 



V X = -U , Vy = V 

The angle tp between vectors v and v (see Fig. |2J) can be found from the relations 

tamp = -^ = — >iuU. (7) 

Let us now consider the situation when the jet propagates along the path OA = vt 
during the time interval t in the FRF, as illustrated in Fig. 121 At the same time the wave 
front from a point-like perturbation (created at point O) reaches a spherical surface with 
radius OB = OC = cd . Two tangent lines AB and AC show boundaries of the Mach 

n 

region [29( with the symmetry axis OA . This region has a conical shape with opening 
angles 8 determined by the expressions (cf. Eq. (|JJ)) 

sm6 = — = - ~c s . 8 
OA v 

Performing inverse transformation from FRF to CMF, it is easy to show that the Mach 
region is modified in two ways. First, it is no longer symmetrical with respect to the jet 
trajectory in the CMF. The insert in Fig. El shows that the boundaries of Mach wave have 
different angles, 9 + ^ with respect to v in this reference frame. One can interpret 
this effect as a consequence of transverse flow which acts like a "wind" deforming the 
Mach cone along the direction OX. On the other hand, the angles of the Mach front with 
respect to the beam (OZ) axis are not changed under the transformation to CMF. We 
conclude, that due to effects of transverse flow, the Mach region in the CMF should have 
a shape of a deformed cone with an elliptic base. 
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FIG. 2: Mach region created by jet moving with velocity v orthogonal to the fluid velocity u . 
Main plot and insert correspond to FRF and CMF, respectively. It is assumed that jet moves 
from O to A in FRF. Dotted circle represents the front of sound wave generated at point O. 

To find the Mach angles 9± it is useful to introduce the angles if)± of the boundary 
lines AB and AC with respect to the OX axis in the FRF (see Fig. EJ. Under the Lorentz 
boost to CMF, ip± are transformed to angles if>± which can be found by using the relations 

cot tp± = 7„ cot i/j± = 7„ tan (9 ± tp) . (9) 

The final expressions for Mach angles 9± = it/ 2 — ip± take the form 

tan 9± = cot ip± ~ 7„ . (10) 

1 =F IsCs^uU 

The last equality gives the approximate formula in the ultrarelativistic case v ~ 1 j^] . 
According to Eq. (jlOj) . 9_ becomes negative for supersonic flows, i.e. at u > c s . For 
such strong flows, the jet trajectory lies outside the Mach region. One can see that the 
difference of the Mach angles 9± in moving and static matter is approximately linear in u . 




FIG. 3: Angles of Mach region created by a jet moving transversely (solid and dashed curves) 
and collinearly (dashed-dotted line) to the fluid velocity u in the CMF. All curves correspond 
to the case c 2 s = 1/3 . Arrow marks the value u = c s . 

Figure 01 shows the numerical values of the Mach angles for an ultrarelativistic jet moving 
through the QGP transversely or collinearly to its flow velocity. The solid and dashed 
curves are calculated by using Eq. (fTU)) with c s — 1/ V3 . We point out a much stronger 
sensitivity of the Mach angles 9± to the transverse flow velocity as compared with the 
collinear flow. 

To discuss possible observable effects, in Fig. 0] we schematically show events with 
different di-jet axes AiBi (i = 1,2,3) with respect to the center of a fireball In the 
2 — 2' event, the away-side jet '2' propagates along the diameter A2B2 , i.e. collinearly 
with respect to the collective flow. In the two other cases, the di-jet axes are oriented 
along the chords, A1B1 and A 3 B 3 , respectively. In such events, the fluid velocity has 
both transverse and collinear components with respect to the jet axis. In Fig. H]we also 
show how the Mach fronts will be deformed in expanding matter. It is easy to see that 
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FIG. 4: Schematic picture of Mach shocks from jets 1,2,3 propagating through the fireball 
matter (shaded circle) created in a central heavy-ion collision. Dotted arrows represent local 
velocities of fireball expansion. Thick downward arrows show associated trigger jets. The Mach 
shock boundaries are shown by solid lines. Short-dashed lines give positions of shock fronts in 
the case of static fireball. 



the radial expansion of the fireball should cause broadening of the sideward peaks in 
the A0-distributions of associated hadrons. Due to the radial expansion, the peaks will 
acquire an additional width of the order of < 9 + — 6- > . Here 9± are local values of 
the Mach angles in individual events. The angular brackets mean averaging over the jet 
trajectory in a given event as well as over all events with different positions of di-jet 
axes. Assuming that particle emission is perpendicular to the surface of Mach cone and 
taking < u >~ 0.4, c s ~ l/y/3 , we estimate the angular spread of emitted hadrons in the 
range 30° — 50° . This is comparable with the half distance between the away-side peaks 
of the A0 distribution observed by the STAR and PHENIX Collaborations fl J E Q 
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On the basis of this analysis we conclude that in individual events the sideward maxima 
should be asymmetric and more narrow than in an ensemble of different events. Due to 
a stronger absorption of particles emitted from the inner part of the shock (events 1, 3 in 
Fig. H|, the two peaks may have different amplitudes. We think that these effects can be 
observed by measuring three-particle correlations (see next section). 

There is one more reason for broadening of the A(/>-distributions which one should 
keep in mind when comparing with experimental data: due to the momentum spread 
of the initial parton distributions, Ap* < 1 GeV, the di-jet system has a nonzero total 
momentum with respect to the global cm. frame. As a consequence, the angle 9* between 
the trigger- and the away-side jet is generally-speaking not equal to tt , as was assumed 
above. Taking typical momenta of initial partons as po > with po > 4 — 6 GeV 
we estimate the angular spread as |7r — ~ Ap*/p < 0.1 . Therefore, the considered 
broadening should be much less than the typical shift of the Mach angles due to the 
collective flow. 



III. INFLUENCE OF LONGITUDINAL EXPANSION OF QGP 

In the preceding section, only effects of transverse flow have been considered. On the 
other hand, experimental data on the jet-induced pseudorapidity correlations in central 
Au+Au collisions at RHIC energies have appeared recently 25[. According to these data, 
"near-side" hadrons (emitted in the forward hemisphere, A<p < tt/2, with respect to 
the trigger jet) exhibit much broader distributions in the relative pseudorapidity, Arj , as 
compared with pp and d+Au interactions. The widths of the Arj distributions increase 
with decreasing p? of secondary hadrons, reaching values of At] ~ 1 for px — * . 

This broadening can be naturally explained as a consequence of the longitudinal ex- 
pansion of fireball matter created at early stages of a heavy-ion collision. Let us consider 
a trigger jet emitted at the initial stage of the reaction (t ~ 0) with the cm. pseudora- 
pidity 7] ~ , i.e. in the transverse plane z ~ , where z is the beam axis. Below we 
consider the Bjorken-type scenario of nuclear collisions |2J] and assuming that a cylin- 
drical, longitudinally expanding volume of QGP with radius R is formed at proper time 
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t = \/t 2 — z 2 = t , where Tq ~ 0.4 — 0.6fm/c. In the following discussion we disregard 
the radial expansion of the fireball, assuming that the near-side correlations are formed 
during a comparatively short time, r — tq < R/c s , which is needed for the trigger jet to 
leave the fireball. 

Weak, sound-like, perturbations of fireball matter propagate in z-direction with ve- 
locities (u ± c s )/(l ± uc s ) , where u and c s are, respectively, the collective longitudinal 
velocity and the sound speed at a given space-time point (t, z) . To find the trajectory of a 



sound wave front z = z + {t) , one should solve the equation for the C + characteristics 

dz + u + c. 



221: 



dt 1 + uc s 



(11) 

z=z + (t) 



Assuming the scaling law, u = z/t = tanh^, for the longitudinal expansion of the fire- 
ball, one can solve Eq. (fTTj) analytically, even in the case when c s is a function of r . 
Indeed, changing the variables from (t, z) to (r, 77) , one gets the following equation for 
the pseudorapidity r/ + of the sound front 

= c s (r) . (12) 

dr 

Assuming that the sound wave is excited at t = r , z = , we get 

V+(r) = J ^c s (t). (13) 

TO 

The expression for the C_ characteristics, 77 = t)-(t) , is given by the r.h.s. of (|T3*jl with an 
additional negative sign. From this derivation we conclude that at fixed proper time t > tq 
a primary jet, created at r ~ 0, r\ = , disturbs the fireball matter in the pseudorapidity 
region r/_ < rj < r] + . In the case c s = const one obtains a simple expression 

7]± = ±c s \og (?-J . (14) 

Within this approach one may expect nonzero correlations of secondaries within the in- 
terval of relative pseudorapidities r/_ ^ Ai] ^ i] + . Substituting r = 0.5 fm/c, r = 5 fm/c, 
c s = l/-\/3 we get an estimate |A?]| max ~ 1.3. This agrees quite well with the experi- 

B 



mental data of Refs. |2j, |25j (after subtraction of the "background" pp peak from the 
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experimental A<p — Arj distribution). As one can see from Eq. (|13j). the width of the Arj 
distribution is sensitive to the QGP formation time To and, therefore, measurements of 
these distributions in different pt intervals can be used as a clock for the thermalization 
process. More detailed information about the Arj — A<p asymmetry can be extracted from 
3-particle correlations. In particular, one can use such correlations to select events where 
the trigger and two associated hadrons have momenta in the axial plane (containing the 
beam axis). If the pseudorapidity of the trigger jet is small, both associated hadrons will 
probe the longitudinal flow of the QGP. 

In the case of a strong shock we expect that the matter in the inner part of the Mach 



region will be swept to its surface 



32]. This will result in a local depletion of parton density 



inside the Mach cone. Such a behavior should be observed as a dip in pseudorapidity 
density of the associated hadrons around r] = . It seems that the experimental data on 
near-side production Q indeed show such a dip for soft secondaries. 



IV. CONCLUSION 



In this paper we have investigated properties of Mach shock waves induced by high- 
energy partons propagating through dense quark-gluon matter created in heavy-ion colli- 
sions. By using simple kinematic and hydrodynamic relations we have analyzed deforma- 
tions of these shocks due to the radial and longitudinal expansion of QGP. Taking typical 
flow parameters expected in central collisions of nuclei at RHIC and LHC energies we 
show that the shape and orientation of Mach regions are strongly modified as compared 
to the case of static (nonexpanding) medium. This may obscure observable signatures of 
Mach collective waves. 

In the future we are going to take into account several additional effects, in particular, 
bending of the leading parton's trajectory due to the friction force from moving medium. 
In principle, this bending may also produce isolated peaks in the away-side A0 distribu- 
tions. On the contrary, the Mach waves should generate ring-like maxima around the jet 
axis. We believe that measurements of three-particle correlations with a broad coverage 
of angular space can be used to differentiate between these two possibilities. We are also 
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planning to investigate properties of strong Mach shocks, in particular, their sensitivity 
to the equation of state of QGP and to effects of collective expansion. 
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